X-RAY AND RADIO EMISSION FROM 
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Abstract 

The collision of the hypersonic winds in early-type bi- 
naries produces shock heated gas, which radiates ther- 
mal X-ray emission, and relativistic electrons, which 
emit nonthermal radio emission. We review our cur- 
rent understanding of the emission in these spectral re- 
gions and discuss models which have been developed 
for the interpretation of this emission. Physical pro- 
cesses which affect the resulting emission are reviewed 
and ideas for the future are noted. 

1 Introduction 

It is well established that massive O-type and Wolf- 
Rayet (WR) stars have powerful, highly supersonic 
winds with high mass-loss rates and terminal speeds. 
These stars are all hot and highly luminous for their 
mass and their winds are driven by momentum trans- 
ferred from the radiation field by line scattering. Rep- 
resentative values for some of the defining parameters 
of massive stars are noted in Tabled 

In addition, a high proportion of massive stars are in 
binary and multiple systems; a volume-limited sam- 
ple found more than 40% of WR stars in multiples 
( van der HuchtL 2001 ). A variety of phenomena re- 
lated to the collision of the two winds can be ob- 
served. For instance, infrared (IR), optical and ultravi- 
olet spectra of many WR binaries show phase-locked 
varia bility in P-Cygni and flat-topped emission lines 
(e.g., Wises & Gies, 1993; Rauw. Vreux & Bohanna: 
[l99l: IStevens & Howatthl 11994 iHilletall I2OO1 
these are signatures of wind-wind collisions. Dust 
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formation may also occur in the wind collision re- 
gion (WCR) and is observed as an IR excess. Bi- 
nary systems which exhibit dust formation are catego- 
rized either in to persistent or episodic dust makers (see 
Willi amsLl2002L for a review). 

Excess X-ray emission from shock-heated plasma pro- 
duced in the WCR is another wind signature - in 
some systems the X-ray luminosity is two orders of 
magnitude higher than suggested by the canonical re- 
lationship, L^/Lh n \ ~ 10~ 7 , from early-type star s 
( Chlebowski & GarmanvL Il99ll: iMoffat et all l2002lS . 
Early-type binaries are also found to be stronger X- 
ray sources than the i r single counterparts in a statis- 
tical sense JPollockL Il987l: iQilebowski & Garmany, 
1991). Massive binaries often display phase-locked 
orbital X-ray variability res ulting from the changing 
line of sight into the sy stem ( Willis. S child & Stevensl 
1 19951 Icorcoranl Il996). Recent Chandra observations 
indicate that the X-ray emission is nonisothermal, that 
the lines are in general narrow, and that there is strong 



2001 



emission from forb idden lines ( Cor coran et al 
Poll ock et aill2004h . 

Finally, nonthermal radio emission from WR stars is 
also a good indicator of win d-wind interaction (e.g., 
Dougherty & Williams! 2000h . In binary systems the 



strong shocks formed by the collision of the two stel- 
lar winds are natural sites for particle acceleration. 
Spatially resolved emission from the WCR in sev- 
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Williams etall 1 19971: iNiemela et all Il998h . has con- 
firmed this picture. 

In this paper we focus on the theoretical modeling 
of X-ray and radio emission from colliding wind bi- 
naries (CWB). X-ray observations provide a direct 
probe of the conditions within the wind-wind colli- 
sion zone and of t he unshocked attenuating material 
in the system (e.g ., Stevens . Blondin & PollockLll992l : 
Pittard & Stevens!. 1997h . Analysis of the X-rav emis- 
sion can determine the mass-loss rates of the stars and 
the speeds of the stellar winds. It is also possible to 
quantify various processes which occur in massive stel- 
lar binaries and high Mach number shocks. In contrast, 
the analysis of nonthermal radio emission from CWB 
systems can yield estimates of the efficiency of parti- 
cle acceleration and the equipartition magnetic fields 
at the surface of the stars, and the mass-loss rates and 
clumping fractions of the stellar winds. CWB sys- 
tems are also excellent objects for investigating the 
physics of diffusive shock acceleration as their geom- 
etry is in many cases very simple, and because they 
provide access to a higher density regime than is pos- 
sible through studies of supernova remnants. Systems 
with highly eccentric orbits which show large mod- 
ulation of the nonthermal emission, such as WR 140 

1 Williams et ailll99(ll:lwhite & BeckerMl995h . should 
be especially useful in this regard. 

This paper is organized as follows. In Section |2] we 
briefly review the basic geometry of the wind-wind 
interaction. In Section PJ] we review the current state 
of the theoretical models used to analyze the X-ray 
emission from colliding wind systems. Systems where 
the cooling lengths may be resolved are amenable to 
study with hydrodynamical codes, and have tradition- 
ally been the focus of work in this area. However, a 
recent development in the modeling of highly radia- 
tive systems is reported. Models which focus on the 
radio emission are reviewed in Section |4l with partic- 
ular attention paid to the first hydrodynamical-based 
investigation which includes thermal and nonthermal 
processes. Section |3] contains a brief review of other 
mechanisms related to CWBs, and in Section |6] we 
summarize and conclude. 

2 Geometry of the wind- wind collision 

Consider the wind-wind collision in a wide binary as 
shown in the schematic in Fig. ^ The winds collide 
at high speeds (v ~ 1000 — 3000 km s^ 1 ) creating a 
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Figure 1: Schematic figure of the geometry of the wind- 
wind collision in a wide binary. For the purposes of this 
work we assume that the binary consists of a WR and OB- 
star, though other combinations (such as O+O-star) also oc- 
cur. Since the binary is assumed to be wide the winds collide 
after they have reached their terminal speeds, and the distor- 
tion of the collision region due to the orbital motion of the 
stars is not visible until far downstream. The half-opening 
angle is indicated. Note that the stars are enlarged for clar- 
ity. 



region of high temperature shocked gas between them 
(T ~ 10 7 - 10 8 K). The position of the WCR is deter- 
mined by ram-pressure balance requirements, and the 
shocked winds are separated by a contact discontinuity. 
If the winds are of equal strength the contact disconti- 
nuity will be a plane midway between the stars, but 
in the more general case the WCR wraps around the 
star with the weaker wind. Maximum values for the 
density, pressure, and temperature occur at the stag- 
nation point, a mathematical singularity which occurs 
where the contact discontinuity intercepts the line be- 
tween the centers of the stars. If the winds are spherical 
the distances r\yR and tob from the WR and OB-stars, 
respectively, to the stagnation point are 



1 



rwR 



1+rJ l/2 D > 1+?? l/2 

where the wind momentum ratio, 

Mob^ob 



1/2 



V 



D, (1) 



(2) 



D is the stellar separation, Mwr and Mob are the 
mass-loss rates of the WR and OB stars, and u\vr and 
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Table 1: Parameters of different evolutionary stages during the life of massive stars. Very massive stars (M > 60 M Q ) are 
thought to evolve from an O-type main sequence star, through a luminous blue variable (LBV) phase to a Wolf-Rayet (WR) 
star, before ending their lives in a supernova (or hypernova) explosion. These stages represent the transition from core-H 
burning to shell-H burning to core-He burning. 



Phase 


Lifetime (yrs) 


MtMoyr- 1 ) 


Woo (km 


s" 1 ) 





~few x 10 B 


~ 10" y - 1(T 5 


~ 1000 


- 3000 


LBV 


~ 10 4 


~ 10~ 3 (+outbursts) 


~ 300- 


500 


WR 


~few x 10 5 


~ 10" 5 - 10" 4 


~ 1000 


-4000 



vob are their wind speeds. From the values in Table 12 
we see that the dimensionless parameter 77 is usually 
small i.e. the WCR is nearer to the OB star than to the 
WR star. 

This simple picture is a good approximation in wide 
binaries, where the winds are accelerated to their ter- 
minal speeds before they collide. The half opening an- 
gle is the angle between the line-of-centers through the 
stars and the asymptotic direction of the contact dis- 
continuity (in the absence of orbital motion). In bina- 
ries where the winds achieve terminal speeds the half- 
opening angle of the WCR (specifically the contact dis- 
continuity) approaches an asymptotic value which is 
approximated by 
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(3) 



for 10 -4 < 77 < 1 (see lEichler & UsovLll993l) . The an- 
gle between the contact discontinuity and each shock 
is largest when the hot gas cools primarily through adi- 
abatic expansion as it flows out of the system (here- 
after "adiabatic systems"), and approaches zero if the 
WCR is able to cool efficiently (hereafter "radiative 
systems"). 

Whether the system is radiative or adiabatic can be de- 
termined by the characteristic coo ling parameter (see 
Stevens. Blondin & Pollock! Il992h 



X 



^cool 



yu 



M_ 7 ' 



(4) 



where i coo i and t^ yn are characteristic timescales for 
cooling and for the flow dynamics, u§ is the wind speed 
in units of 1000 km s -1 , D\2 is the stellar separation 
in units of 10 12 cm, and M_7 is the mass-loss rate of 
the star in units of 10 -7 M Q yr" 1 . Here we have as- 
sumed that the temperature of the shocked gas corre- 
spond roughly to the minimum in the cooling function 
at T ~ 10 7 K (see Fig.|2]for sample cooling functions). 



Values of x and other interesting parameters for some 
of the most well-studied systems are listed in Table |2] 
It is clear from Table |2] that CWBs are a very diverse 
class of objects, and this is reflected in the wide range 
of phenomena that occur in them. 

3 X-ray emission from CWBs 



Prilutskii & Usov J 19761) and Icherenas hchiikl Jl976) 
proposed that binary star systems, in which one or both 
components are massive OB or WR stars with strong 
stellar winds, should emit a detectable flux of X-rays 
as a result of the shocked gas formed where the winds 
collide. 

We see from Eq.|4]that since x oc Z), short-period sys- 
tems tend to be radiative, while wide systems tend to 
be adiabatic. In highly radiative systems (x <C 1), the 
X-ray luminosity L x oc Mv 2 . If both of the shocked 
winds satisfy x <C 1 then the faster wind w ill dom- 
inate the emission (Pittard & S tevensL 2002h - other- 
wise, the wind with the lower value of x will i n gen- 
eral dominate L x . Alternatively, if both of the shocked 

winds are largely adiabatic, then L x oc M /v 3 ' 2 D 
dStevens. Blondin & Pollock! 1992), and the stronger 



wind dominates the emission. In this case, > yuyo 
of the X-ray luminosity is emitted within an off -axis 
distance equal to D (see Pittard & StevensL 2002 , and 
references therein). 

The particular method used to construct a theoretical 
model of the X-ray emission from a CWB is to a large 
extent determined by the nature of the WCR. Most 
work to date has focused on systems where the stellar 
separation is fairly large and where the WCR is largely 
adiabatic. In such systems it is normal to use a hydro- 
dynamical model to obtain a description of the WCR 
which can then be used to calculate synthetic X-ray 
data. We review the steps behind the construction of 
these models in Sec. 13. II 
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Table 2: Parameters of interest for some of the most important colliding wind systems. The systems are ordered 
according to their orbital period. Where the orbit is eccentric a range in the possible values is noted. The WCR 
may crash onto the O-star photosphere in V444 Cyg, which makes an estimate for xob somewhat difficult in this 
system. 



System 


Orbital period (d) 


Separation (AU) 


Density ( cm 3 ) 


XWR 




WR 139 (V444 Cyg) 


4.2 


0.2 


~10 1U 


< 1 


? 


WR 11 ( 7 2 Vel) 


78.5 


0.8-1.6 


~10 9 


~ 0.5 - 1 


250 - 500 


WR 140 


2899 


~ 1.9-30.7 


~ 10 9 - 10 7 


~ 2 - 50 


100 - 1000 


WR 147 


> 10 5 


> 410 


< 10 4 


> 30 


> 1000 



In short-period systems the WCR is often highly ra- 
diative. Due to the difficulties involved in construct- 
ing hydrodynamical models under these conditions, to 
date there have been few attempts to calculate the X- 
ray emission from such systems. In Sec. l3.2l we review 
recent work which takes a distinctly different approach 
which bypasses problems which a traditional hydrody- 
namical simulation encounters in this regime. 

3.1 Hydrodynamical models of the WCR 



were presented 



Early numerical models 
by (a mongst others) Lebedev & Mvas nikov 
(198 81). Euo. McCrav & Mac Low! Jl990h . and 
Steve ns. Blo ndin & Pollockl(ll992h . The main ingredi- 
ents involved in these models are summarized in Fig.|2] 
A hydrodynamical code is used to model the dynamics 
of the WCR. In this example the deflection of the 
WCR by the orbital motion of the stars is assumed to 
be negligible, the WCR remains axis-symmetric, and a 
2D (r, z) calculation is performed. A small amount of 
artificial viscosity is used in the calculations to prevent 
the development of the odd-even decoupling (i.e. 
carbuncle) instability at the apex of the collision zone. 
This numerical artifact develops when sh ocks are 
locally aligned with th e hydrodynamic grid dWalden. 
19941 : Le Veauel Il998h . In Fig. E rj = 0.1, and the 
stellar winds were assumed to have reached their 
terminal speeds. We also specified uwr = vo, so 
the velocity shear along the contact discontinuity 
is small and there is no noticeable development of 
the Kelvin-Helmholtz instability (instabilities are 
discussed in greater detail in Sec. 15.31 . In systems 
where the WCR is largely adiabatic, it is not strictly 
necessary to compute the rate of energy loss of the 
shocked gas through radiative cooling (since by 
definition this is small). However, it is often included 
in order that models are self-consistent. The hot 



plasma is generally assumed to be optically thin and 
in collisional ionization equilibrium, in which case 
the cooling rate as a function of temperature may be 
calculated usi ng, for example, the MEKAL p lasma 
emission code (Mewe, Kaastra & Liedahl, 1995). 

Cooling curves for solar abundances (X[H]= 0.705, 
y[He]= 0.205, Z[everything else]= 0.020) and the 
abundances found for the WC8 star in 7 2 Vel (X = 
0.0, Y = 0.628, Z = 0.372) are shown in the 
top right panel of Fig. |2] Line emission dominates 
the cooling at temperatures below 10 7 K. At higher 
temperatures cooling occurs primarily through thermal 
bremsstrahlung (A(T) oc T 1 / 2 ). The enhanced metal 
abundance in the winds of WR stars leads to more ef- 
ficient cooling than is obtained for material of solar 
abundance. 

To obtain synthetic X-ray data the hydrodynamical 
model is processed through a ray-tracing code which 
solves the radiative transfer equation along lines of 
sight through the hydrodynamical grid. Emissivity and 
opacity data for solar abundance material are shown 
in the bottom left and right panels of Fig. 13 respec- 
tively. The emissivity data was again calculated using 
the MEKAL code, and consists of line emission plus 
thermal bremsstrahlung continuum. At T = 10 8 K 
most elements are completely ionized and there is lit- 
tle line emission in comparison to lower temperatures, 
while the bremsstrahlung emission extends to higher 
energies. The opacity data in the bottom right panel 
of Fig. |2] w as calculated u sing the Cloudy photoion- 
ization code dFerlandl 2001 ). Distinct absorption edges 
are clearly visible. As the temperature increases and 
elements become more and more ionized the absorp- 
tion cross-section decreases. 

Once synthetic spectra and images have been calcu- 
lated they can be compared against actual X-ray data, 
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Figure 2: The main ingredients involved in the construction of models of the X-ray emission from colliding wind systems 
which are largely adiabatic. Top left: a density plot from a hydrodynamical model of the wind-wind collision. Top right: 
cooling curves for hot, thermal plasma. Bottom left: emissivity data (for solar abundances). Bottom right: opacity data (for 
solar abundances). 



and the properties of CWBs determined. The most 
rigorous method involves the calculation of a grid of 
synthetic X-ray data from a large number of hydro- 
dynamical models in which the stellar wind parame- 
ters are systematically varied. An XSPEC table model 
of the synthetic X-ray data grid can then be gener- 
ated, and fitted against actual observational data. Con- 
straints on the stellar mass-loss rates and wind terminal 
speeds have b e en obt ained by this method for 7 2 Vel 
( Stevens et all Il996h . WR 140 dzhekov & Skinner . 
2000h . and r, Tcarinae dPittard & Corcorarl l2002h . In 
each of these works the X-ray derived mass-loss rate 
for the star with the stronger wind is substantially 
lower than the best estimate from radio observations 
at that time. This is largely due to the fact that the X- 



ray derivation appears to be insensitive to wind clump- 
ing (unlike radio-based derivations of M). The X-ray 
based method also enjoys further advantages. For in- 
stance, in the case of r\ Carinae the high energy X-ray's 
are much less affected by the complex circumstellar 
environment than emission at longer wavelengths. An 
additional benefit over other techniques is the ability to 
place constraints on the properties of the star with the 
weaker wind. 

To illustrate this method at work, Fig. |3] shows a fit to 
a Ch andra grating spectrum of r\ Carinae (reproduced 
from Pitta rd & Corcorarl 2002). r\ Carinae is arguably 
the most infamous massive star in our Galaxy, being 
best known as the survivor of the greatest non-terminal 
stellar explosion ever recorded. The ejecta have since 
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formed the beautifu l bipolar nebula kno wn as the Ho- 
munculus (see, e. g.JMorse et allll998l) . Study of this 
nebula (see, e.g., Sm ith et all 120031) plays a central 
role in the quest for a deeper understanding of this 
object. Unfortunately, it seems that astronomers can- 
not have their cake and eat it, since the nebula also 
acts as a screen which significantly obscures the cen- 
tral star. For a long time 77 Carinae was thought to 
be a single massive star and was classified as an lu- 
minous blue variable (LBV), but there is now over- 
whelming evidence for binarit y and a strong wind- 
wind collision in thi s system ( Damine li et all 2000; 
Corcoran et all 120011) . Tantalizing evidence also exists 
to relate the "Great Eruption" in 1843 with the close 
approach of the stars during their passage through pe- 
riastron. Whether the tidal influence of the companion 
star was in part responsible for triggering the enormous 
outburst is one of the most important questions con- 
cerning 7] Carinae today. 

77 Carinae is in fact a highly unusual example of 
a CWB. It may be the only one known in which 
one of the components is an L BV star (although 
HD 5980 may be another - see, e.g. lNaze et all 12002). 
and as such the low wind speed of the primary star 
(v ~ 500 km s -1 ) means that the hard X-ray emis- 
sion (above 2 keV) is generated almost exclusively by 
shocked material from the secondary wind. At the time 
of the Chandra observation (shown in Fig. |3jl the at- 
tenuation of this hard emission by the stellar winds is 
small and the line of sight into the system is through 
the less dense wind of the companion star. This gives 
rise to the unusual situation that the hard X-ray emis- 
sion is not directly sensitive to the mass-loss rate of the 
primary, but only indirectly through the wind momen- 
tum ratio 77. Fortunately, the duration of the lightcurve 
minimum can be used to estimate 77. The best fit to 
the X-ray data then gives M 10 _5 M Q yr _1 for 
the secondary star and v « 3000 km s -1 for its wind. 
With ■>] = 0.2 and assuming v ~ 600 km s -1 for 
the wind of the primary star, the fit in Fig. |3] implies 
M 2.5 x 10 -4 M yr _1 for the primary star. As 
noted earlier, this value is substantially smaller than 
commonly inferred (in the literature there is a large 
scatter in the derived values, though the general view is 
that M ~ 10 -3 M & yr _1 ). However, some problems 
with a mass-loss rate as high as 10~ 3 M Q yr _1 have 
been identified: fits to HST data show absorption com- 
ponents which are too strong an d electron-scattering 
wings which were overestimated (Hillie r et all 12001). 



data and folded 




Figure 3: The best fit to a Chandra grat ing spectrum of 
77 Carinae (from lPittard & Corcoran! [20021) . If we take the 
terminal speed of the primary wind as v\ = 600 kms -1 , 
then the fit implies that M 1 = 2.5 x 10" 4 M yr -1 . 



Likewise, the X-ray derived results face their own set 
of problems. For instance, the derived wind param- 
eters of the secondary star suggest that it is either a 
very massive O-star or perhaps a WR star, but no ev- 
idence for a companion has been found to date in op- 
tical spectra. It will be interesting, ther efore, to see if 
the results in lPittard & Corcoran (2002) are consistent 
with those from further X-ray observations. Analysis 
of additional Chandra spectra taken through periastron 
passage is ongoing, and we hope to report results soon. 

Since the launch of Chandra and XMM-Newton, it has 
also become possible to measure the shift and width of 
individual X-ray lines from CWBs. The high spectral 
resolution offered by the gratings on these telescopes 
allows access to a wide range of important diagnostic 
data which can be used to place additional constraints 
on the wind-wind collision. The theoretical study of 
line profiles in CWBs is in its infancy, although some 
unexpected results have already been obtained. We re- 
fer the interested reader to Henley et al. (these pro- 
ceedings). 

3.2 Modeling X-ray emission from highly radia- 
tive systems 

It is impossible to use the techniques described in the 
previous section when the shocked gas is highly ra- 
diative because in order to calculate the emission the 
cooling length must be spatially resolved. Hydrody- 
namical models (including those with adaptive-mesh- 
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refinement) impose a limit to the spatial resolution 
which can be much larger than the cooling length of 
the gas in such systems. Another problem is that 
the rapid cooling of hot gas creates cold dense sheets 
which are unstable to a varie t y of in stabilities (see, e.g., 
Stevens. Blondin & PollockL 1992b . and which intro- 
duce a large time dependence into the computed emis- 
sion. 

These short-comings have been rec ently addressed by 
Antok hin. Owocki & Brownl (2004). Their solution is 
to separate the calculation of the small-scale (local) 
cooling structure of the post-shock gas from the global 
structure of the wind ram-pressure balance surface, and 
to assume that these structures are time-independent. 
The steady-state shape of the contact discontinuity in 
a CWB has been comprehe nsively discussed in the lit- 
eratur e (see references in lAntokhin. Owocki & Brown , 
2004), though in these works the winds were assumed 
to have reached their terminal speeds before colliding. 
However, as the post-shock gas becomes increasingly 
radiative with decreasing stellar separation, such mod- 
els are most applicable to short-period systems where 
the wind speeds will often have not reached their termi- 
nal values at the interaction f r ont. T o account for this, 
Antokhin. Owocki & Brow n (2004) assumed "beta" 1 
velocity laws in their model. 

In the highly radiative limit, the timescale for shock 
heated material to cool is small compared to the flow 
time for material to advect along a substantial arc of 
the contact discontinuity, and the post-shock gas can 
be assumed to occupy a layer which is geometrically 
thin with respect to other scales in the system (e.g., the 
binary separation, and the radius of curvature of the 
contact discontinuity). The expansion of the shocked 
gas as it flows out of the system normally couples the 
internal evolution of the post-shock gas to the global 
wind-wind structure, since adiabatic cooling can then 
take place. However, with the above assumptions one 
can adopt an "on-the-spot" treatment of the radiative 
emission at each point along the interaction front, and 
the post-shock structure can be simplified to planar ge- 



Contact surface 



'A beta velocity law has the form: 

v(r) =v oc (l-R/rf 

where R is the radius of the star. /3 = 1 is often used in analytical 
works, [3 = 0.8 is fairly representative of O-stars, while WR stars 
have more slowly accelerating winds (which means higher values 
of beta, though such a formulation is probably very poor for WR 
stars). 



Po 




Po. 



o.o 



Figure 4: A schematic of the post-shock density and tem- 
perature str ucture in a steady, planar cooling la yer (repro- 
duced from Antok hin. Owocki & Browi] ll2004l) . with per- 
mission). The shock front occurs at Iq. 



ometry. 

To obtain the structure of the cooling layer, 

Antokhin. Owocki & Brown (2004) adopted an iso- 
baric approximation. Assuming that the energy loss 
in the shocked gas is solely from radiative cooling, the 
rate of temperature decrease in the post-shock gas is 
given by, 



5k dT A trns 

-pv— = -n e n H A(Tj, 



2fim p dl 



(5) 



where / is a distance from the shock front, and A(T) is 
the cooling function. Equation [5] also gives the width 
of the cooling layer, and a schematic of the post-shock 
structure is shown in Fig. |4] 

An example of the global structure of the WCR, to- 
gether with the width of the post-shock cooling layer 
either side of the contact discontinuity, is shown in 
Fig. 13 An initial application of this method has been 
made to an XMM-Newton observation of the short- 



period binary HD 159176 dDe Becker et all |2004), 



with good success achieved in fitting the spectral 
shape. However, in this work it was not possible to ob- 
tain a good fit which was simultaneously in agreement 
with the predicted luminosity from the model. 

This method shows great promise and it will clearly 
be useful to develop it further. Future work on this 
model could include extension to three dimentions 
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Figure 5 : Th e structure of the wind-wind intera ction from 
the model of lAntokhin. Owocki & Browi] J2004I) (repro- 
duced with permission). The solid line marks the con- 
tact discontinuity, while the dashed lines indicate the shock 
fronts. 



(which would allow deflection of the interaction re- 
gion into a spiral form as a result of orbital motion 
of the stars), and the dyn amical role of the stellar 
radiation fields (see, e.g., Stevens & Pollocil ,1994; 
iGaylev. Owocki & Cranmei , ll997l and further discus- 
sion in Sec. [5J. Further tests against a variety of short- 
period CWBs also need to be made, and in this re- 
spect XMM-Newton data from the archetype system 
V444 Cygni scheduled for May-June 2004 are ideally 
suited. It should also prove illuminating to compare 
the results from this method and those from hydrody- 
namical models in the transition region where \ ~ 1- 



4 Radio emission from early-type stars 

Radio observations of early-type stars reveal that 
they can emit both thermal and nonthermal radia- 
tion. The former is readily explained as arising 
from free-free emi s sion w ithin the strong stellar wind 
( Wright & BarlowL Il975l) . The winds of early-type 
stars are partially optically thick at radio wavelengths, 
and the optical depth at radius r is given by 



r u (r) 



10" 




9u 



1 



m 2 T3 /2 U 2 r 4 



dr.. 



(6) 

where Z is the rms charge of the atoms (Z = 1 for 
singly ionized gas), j e is the number of free electrons 



per ion, ^ is the mean atomic mass of the ions in 
units of mu, g u is the gaunt fac tor and T is the tem- 
peratu re of the wind (see, e.g., Lamers & Cassinelli, 
1999t) . The high mass-loss rate of massive stars means 



that the photospheric radius at radio wavelengths is al- 
ways much larger than the radius of the underlying star. 
For an isothermal O-star wind with T = 10, 000 K, 
M = 10" 6 Moyr" 1 , = 2000 kms" 1 , and con- 
sisting mainly of ionized hydrogen, the r„ = 1 sur- 
face occurs at radii of 1350, 1900 and 4200 R at radio 
wavelengths of 3.6, 6 and 20 cm respectively. For WR 
stars, these radii are even larger. Since the wind speed 
is assumed to have reached its terminal value at these 
distances, radio observations have proved to be a use- 
ful tool for measurements of mass-loss rates (although 
as discussed in Sec. B.ll these suffer from uncertainties 
due to the clumpy nature of such winds). 

The nonthermal emission, on the other hand, is at- 
tributed to synchrotron radiation from relativistic elec- 
trons moving in a magnetic field. The consensus is that 
the relativistic electrons in the winds of massive stars 
are created through diffusive shock acceleration. In 
this process the electrons gain energy through the first- 
order Fermi mechanism as they cross and re-cross the 
shock unt il they are e ventually scattered downstream 
(see, e.g., Bell 1978h . Electrons with Lorent z factor 



7 < 10 5 are frozen into the post-shock flow {white, 



1985). The resulting energy distribution of the rela- 



tivistic electrons has a power law form, i.e. the number 
of relativistic elections per unit volume with energy be- 
tween 7 and 7 + d'y, N(j) oc 7~ n . The energy index, 
n, depends on the shock compression ratio, \c, through 



1223) 



n 



Xc + 2 



Xc 



1 



(7) 



In the case where the shocked plasma radiates effi- 
ciently, the shock is quasi-isothermal, and the com- 
pression ratio can become very large, leading to n ~ 1. 
If instead the shock is adiabatic, then Xc = 4 (for a 
strong shock), and n = 2. 

As already noted, hydrodynamical shocks are a fea- 
ture of the winds of both single massive stars and their 
binary counterparts. In the case of single stars, the in- 
trinsic instability of the radiatively driven wind leads to 
the formation of s hocks which are distr ibuted through- 
out the wind (see Ru nacres & Owocki 2002L and ref- 
erences therein). While in principle these can accel- 
erate electrons to high energies, the intense stellar ra- 
diation field causes significant energy losses through 
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inverse Compton (IC) scattering and electrons accel- 
erated in the inner regions of the stellar wind cannot 
survive out to large radii. Those electrons which are 
responsible for the nonthermal radio emission must 
therefore be accelerate d at large distances from the star 
I Chen & WhiteL 1994). The large photospheric radii 
of massive stars at radio wavelengths imposes a simi- 
lar constraint, since for nonthermal emission to be ob- 
served at radio wavelengths, the relativistic electrons 
must exist outside the radio photospheric radius of the 
stellar wind. 

In binary systems the strong shocks formed by the col- 
lision of the two winds is a natural place for diffu- 
sive shock acceleration of nonthermal electrons. An- 
other possibility concerns the current sheets which re- 
sult from the collision of mag netized stellar winds 
I Jardine. Allen & PorlockLll996h . ' The compression of 
magnetic field lines in the current sheets create en- 
hanced local field strengths, and the resulting electric 
fields may accelerate particles to high energies. Irre- 
spective of which is the dominant acceleration mech- 
anism, wide binary systems have the significant ad- 
vantages that the wind collision may exist well out- 
side the radio photosphere of the stars, and that the 
stellar radiation energy density is substantially di- 
luted in the vicinity of the WCR. This is reflected 
in strong observational support for a CWB interpre- 
tation for all synchrotron em ission from WR stars 
( Dougherty & Williams!. l2000h. though the case is not 
so clear for O-stars (see, e.g.. lRauw et al. , 2002 ). 
However, when synchrotron emission is not detected, 
that does not mean it is not being emitted. Nonthermal 
radio emission generated in the WCR of a short period 
binary may be buried so deep under the radio photo- 
sphere that the heavy opacity prevents its observation. 
For the nonthermal emission to be highly visible we 
require that the binary is wide enough that the apex 
of the WCR is far enough away from the stars that it 
exists outside of their radio photospheres. This separa- 
tion implies a value of the orbital period of order of a 
month. Therefore, there is a selection effect that all the 
CWBs with observed synchrotron are wide binaries. 

4.1 Modeling the radio emission from CWBs 

Most analysis of the radio spectra from colliding wind 
binaries has hitherto been based on a highly simplified 
model consisting of a thermal source plus a point-like 
source of nonthermal emission attenuated by free-free 



absor ption (e.g., Chapman et al. , 19991 : Monnier et al 
2002h . In this model, at frequency u, the observed 
flux (S*° bs ) is related to the thermal flux (Sl h ), the syn- 
chrotron emission arising from the WCR (5^ yn ), and 
the free-free opacity ( rjj) of the circum-system stellar 
wind envelope by 



s: 



obs 



(8) 



From simultaneous measurements in a number of fre- 
quency bands it is possible to solve Eq. [8] to deter- 
mine the relative contribution of thermal and nonther- 
mal emission to the total flux. On the other hand, both 
5^ yn and may display significant time variability. 
The former is expected if the stellar separation varied 
as in an eccentric binary, while the latter is in general 
anticipated since the optical depth to the synchrotron 
emitting region will change as the line of sight sweeps 
around the system. 

While Eq. [8] may at times be adequate, a single- 
valued free-free opacity determined along the 
line-of-sight to a point-like synchrotron emission 
region is a clear over-simplification: observa- 
tions reveal an extended region of synchrotron 
emission 



emiss i on from both WR 146 JDoughertv et al. , 
1994 boughertv. Williams & PollaccoL l2000h and 
WR 147 dMoranetallll989l : IChurc hwell et all ll992l : 
IWilliams et all Il997h - see also Sec. 14.1.11 Further- 
more, use of such over-simplified models has often 
resulted in difficulties in fitting observational data. 
Nowhere are these problems better illustrated than 
for the wonderful dataset which exists for WR 140, 



the a r chetype wide CWB (se e, e.g., IWilliams et al 
1990; White & Becker, 1995). We note that the 



effect of binarity on the thermal radio emiss ion from 
massi ve stars has also been investigated (iStevensl. 
1995b . though this work did not examine nonthermal 
emission. 

Complex magnetohydrodynamic models which in- 
clude the transport of nonthermal electrons have 
been computed to predict the synchrotron emission 
from other astrophysical obje cts, such as radio jets 
( Tregilhs. Jones & RyuL l2004h . However, less atten- 
tion has been focused on models for CWBs, and the re- 
sulting models are less well developed. In fact, CWBs 
are nice laboratories for particle acceleration because 
the fundamental flow parameters of the winds are well 
known. 

As a first step towards the construction of more realis- 
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Frequency (GHz) 



Figure 6: Spectra from the standard model with 0° inclina- 
tion - free-free flux (dashed), magneto-bremsstrahlung flux 
(dotted), and total flux (solid). IC cooling, coulombic cool- 
ing, the Razin effect and synchrotron self absorbtion (SSA) 
are not included in this calculation. 



tic models, Doug herty et al.l (|2p03) obtained a simple 
representation of the spatial distribution of the free-free 
and nonthermal emission from CWB systems through 
the use of a hydrodynamical model of the WCR. In this 
section we note the main results of this work. 

The process of calculating the theoretical radio emis- 
sion from a hydrodynamical model of a CWB is as 
follows. First, the temperature and density values 
on the hydrodynamic grid are used to calculate the 
free-free emission and absorption coefficients from 
each grid cell. The synchrotron emission and syn- 
chrotron self-absorption from each cell within the 
WCR is calculated assuming that the distribution of 
relativistic electrons in each cell can be specified by 
N{^)d^ = C^~ v dq, where p = 2. The stan- 
dard assumption that C is proport ional to the ther- 
mal energy density is made (c.f. IChevaheiL 19821 : 
Miodusze wski. Dwarkadas & Ball 200lb . and it is as- 
sumed that p is spatially invariant over the volume of 
the WCR. The magnetic field in the WCR is assumed 
to be highly tangled, and the synchrotron emissivity 
(P v oc i/ _ (p -1 )/ 2 ) is assumed to be valid over all fre- 
quencies, v. A radiative transfer calculation is then 
computed to obtain the flux and intensity distribution 
at a speci fied frequency. The i nterested reader is re- 
ferred to lDoughertv et all (120031) for further details. 
We shall first discuss the results obtained from a "stan- 
dard" CWB model with the parameters noted in Ta- 
ble [3] These wind values are typical of WR and OB 
stars, and with the adopted binary separation, D, the 



Tabl e 3: Parameters for the "standard" CWB model in 
Dougherty et al. I j2003n . 



Parameter 



Value 



M W R(M yr- 1 ) 2 x 10" 5 

M o (M yr- 1 ) 2 x 10" 6 

«oo,WR (kms^ 1 ) 2000 

vjp (kms" 1 ) 2000 

D(cm) 2 x 10 15 



WCR is largely adiabatic. At a frequency v = 5 GHz, 
D is approximately 10 x the radius of the Tg = 1 sur- 
face of the WR wind. For inclination angles ~0°, the 
lines of sight to the WCR are then optically thin, per- 
mitting investigation of the emission from the wind- 
collision region in the absence of strong free-free ab- 
sorption from the circum-binary stellar wind envelope. 
Solar abundances for the O-star and WC-type abun- 
dances for the WR-star (mass fractions X = 0, Y = 
0.75, Z = 0.25) are adopted in the model. Tempera- 
tures of 10, 000 K and an ionization structure of H + , 
He + and CN0 2+ are assumed for the unshocked stel- 
lar winds, and the system is assumed to be positioned 
at a distance of 1 .0 kpc. Finally, it is assumed that the 
magnetic energy density, Ub, the relativistic electron 
energy density, U ie \, and the thermal energy density, 
U t h, are related by 



B 2 

Ub = — = CB^th, 

8"7T 



(9) 



and 



U 



rel 



n('j)-fmc 2 d'y = Cretan, (10) 



where £b and £ re i are constants, and C/ t h = pri where 
P is the gas pressure and F is the adiabatic index (as- 
sumed to be 5/3, as for an ideal gas). The common as- 
sumption of equipartition corresponds to £b = Crei = 
£. The nonthermal electrons lose energy due to syn- 
chrotron emission and IC scattering in the same ratio 
as the magnetic field energy density to the photon en- 
ergy density, 



P 



sync 



compt 



Ub 
U ph 



(11) 

< T'compt, and 



For all reasonable values of (b, T'sync 
IC cooling is the dominant energy loss mechanism. 
Synchrotron cooling is therefore ignored. 

The maximum post-shock density, temperature, and B- 
fleld occur at the stagnation point between the stars and 
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Figure 7: Intensity distributions from the standard model 
at 0° inclination at 1.6 GHz (top) and 22 GHz (bottom). 
Neither IC cooling, coulombic cooling, the Razin effect or 
SSA, are included in this calculation. The stellar winds and 
the WCR are clearly visible. 



are listed in Table|4] The synchrotron emission arising 
from a single relativistic electron w ith Lorentz factor 7 
is (c.f. Rybicki & Lightmarl [l97^ 



Ply) 



\/3q 3 B since 



m e c" 



F(y/v c 



(12) 



where q is the electron charge, B is the magnetic field 
strength, a is the pitch angle of the particle relative to 



Tabl e 4: Values from the "standard" CWB model in 
Dougherty et al ] <l2003h . 



Parameter 



Value 



"max 


( cm 3 ) 


4 x 10 b 


T 

± max 


(K) 


2 x 10 8 


^max 


(niG) 


6 



the direction of the B-field, m e is the mass of the elec- 
tron, c is the speed of light, and F{u/u c ) is a dimen- 
sionless function describing the total power spectrum 
of the synchrotron emission, with v c the frequency 
where the spectrum cuts off, given by 



3^/ 2 qBsma 
47rm P c 



(13) 



in 



Valu es for F(v/y r ) are tabulated 

Ginzburg & Svrovatskiil (1 19651) and we assume 
that sin a = 1. At the stagnation point z/ c = 0.25 GHz 
(for 7 = 10 2 ) and 25 GHz (for 7 = 10 3 ). 

The spectra and intensity distributions obtained from 
the standard model with ( = 10~ 4 are shown in Fig- [H] 
and Fig.Qrespectively where, for simplicity, only free- 
free absorption is included. Magneto-bremsstrahlung 
emission dominates the total flux below 1 GHz. The 
turnover seen near 400 MHz is due to free-free ab- 
sorption from the unshocked stellar winds. Above 
400 MHz, the magneto-bremsstrahlung spectrum is a 
power-law with a spectral index (S u oc v a ) a = —0.5, 
as expected for a p = 2 energy distribution. The 
thermal spectrum spectral index of +0.6 is also as ex- 
pected, though a slight departure from a power-law ex- 
ists at low frequencies due to the finite size of the hy- 
drodynamical grid. The intensity distributions shown 
in Fig. |7] reveal that at 1.6 GHz the WCR is brighter 
than the stellar winds, while the opposite is true at 
22 GHz. 

Dougherty et al. J2003I) also explored how the radio 
flux varied with stellar separation, and we repeat the 
main conclusions here. We shall first consider how 
varying the stellar separation affects the free-free opac- 
ity to a specific point in the WCR. For a stellar wind 
with an r~ 2 radial density distribution the line-of- 
sight opacity r, at frequency v through the wind is 
r oc cj -3 ^ -21 , where £ is the impact parameter. 
Since £ oc D for a given inclination, the turnover fre- 
quency for a constant opacity value is v oc D~ 10 / 7 . 
Dougherty et al.1 (120031) showed that this relationship 
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is in excellent agreement with the data shown in their 
Fig. 5. 

Now let us consider how the intrinsic synchrotron 
emission depends on the stellar separation. In the 
absence of synchrotron self absorbtion (SSA) or the 
Razin effect, the synchrotron luminosity increases as 
the separation decreases, due to the increased thermal 
energy density in the collision region. In this sce- 
nario, the intrinsic synchrotron emission per unit vol- 
ume P{v) oc £ 3 / 4 n 3 / 4 i/ -1 / 2 for an electron power-law 
index p = 2. Since the post-shock density oc D~ 2 
and the volume of the emitting region scales as D 3 , 
the total synchrotron emission from the entire wind 
collision volume scale s as D~ l l 2 v~ l l 2 , as noted by 



Table 5: Parameters for the WR 147 model in 



Dougher ty et al.l d2003r) and displayed in their Fig. 6. 




Frequency (GHz) 



Figure 8: The radio spectrum of WR 147. The 
deduced synchrotron emission is represented by solid 
data points, while the total flux is represented by the 
hollow data points. Estimates from different authors 
are as follows: solid squares (two separate e stimates 
from the same observational data, Skinner etall Il999h , 
solid pentagons jsetia Gun awan et all l2Q0l (). open circle s 
dChurchwell et all Il992l IContreras & Rodnguezl Il999h . 
The lines are the total (solid), thermal (dotted) and syn- 
chrotron (dot-dashed) spectra of the i = 0° model with 
C = 7.03 x 10~ 3 , where SSA, Razin and free-free absorp- 
tion are included in the radiative transfer calculations. The 
model winds were assumed to be clumpy with a volume- 
filling factor / = 0.134. Models at +30° and -30° are also 
consistent with the data. 



4.1.1 Modeling the radio emission from WR 147 

WR 147 is notable for being among the brightest 
WR stars at radio frequencies, and for being one of 
two systems in which the thermal and synchrotron 
emission are observed to arise from two spatially re- 



Doughert v et all d2003). 



Parameter 


Value 


M W N8 (Moyr- 1 ) 


2 x 10" 5 


Mob (Meyr- 1 ) 


3.8 x 1(T 7 


^oo,WN8 (kms -1 ) 


950 


Woo.OB (kms -1 ) 


1000 



solved regions (see, e.g., Willia ms et al.ll 19971 and ref- 
erences therein). Furthermore, it is one of a hand- 
ful of WR+OB binary systems where the two stars 
are resolved into a v isual pair (iNiemela et al. , 1998; 



Willi ams et all 119971) with a projected stellar separa- 
tion given by Dcos i = 0.6 35 ± 0.020 / . / . At the esti- 
mated distan ce of ^ 0.65 kpc JChurchwelletallll992l : 
iMorris et all 120001) this corresponds to a separation 
D ~ 415/ cos i AU, and the relationship between D 
and i represents an important constraint for any mod- 
els of the system . As th e inclination angle is unknown, 
iDoughertv et al.ld2003h investigated models for several 
different values of i, which in turn requires different 
values of the physical separation, D, to maintain the 
observed angular separation. 

The radio spectrum of WR 147 is perhaps the best ob- 
served of all massive binary systems, with radiome- 
try extending from 353 MHz to 42.9 GHz. However, 
differences in excess of 50% exist in the derived syn- 
chrotron flux at some frequencies (cf. Fig. giving 
rise to uncerta inty in the position of th e low frequency 



turnover (see Doughe rty et all I2003L for further de- 
tails). These differences form the fundamental reason 
for conflicting conclusions about the nature of the un- 
derly ing electron energy s pectrum in the curren t litera- 
ture dSkinner et allll999l : lsetia Gunawan et ailbOOlh . 

Full details of the mode ling of the radio spectr um of 
WR 147 can be found in Doug herty etall d2003l) . Val- 



ues for the mass-loss rates of the stars, and the termi- 
nal velocities of the winds are noted in Table |5] These 
were taken from the literature, as were details of the 
wind compositions. The temperature for both stellar 
winds was assumed to be 10 kK, and the dominant ion- 
ization states were taken to be H + , He + and CN0 2+ . 
The thermal flux is only weakly dependent on the as- 
sumed wind temperature (through the Gaunt factor), if 
the ionization state is fixed. 

Models with i = 0, 30 and 60° were computed, and 
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values of ( were estimated for each model by match- 
ing the mean of the "observed" synchrotron 4.86 GHz 
fluxes, which was taken to be 14.1 ± 0.3 mJy. The 
spectrum between 0.2 and 80 GHz was then calculated 
assuming this constant value of £. To fit the thermal 
flux with the chosen value of Mwns the winds were 
required to be clumped. 

In Fig. |8] the resulting total, synchrotron and thermal 
spectra from one of the models (i = 0, and includ- 
ing SSA, the Razin effect and free-free absorption) is 
shown against the observed spectra. Though this is not 
the best fit to the data in a formal sense, the model fits 
the data reasonably well given the approximations and 
assumptions which it contains. The shortcomings are 
that the total flux in the 5 to 8 GHz range is under- 
estimated, and the synchrotron emission does not turn 
down somewhere around 10 to 20 GHz as suggested 
by the observed data. The shortfall in the 5-8 GHz 
total flux could be eliminated by increasing the ther- 
mal emission by a few mJy. This could be achieved by 
simply increasing M / in the stellar winds. While 
this would require that the synchrotron emission at fre- 
quencies higher than ~ 10 GHz be lower than in the 
current models, this may well be the case if IC cooling 
is taken into account, as explained below. 

The fit shown in Fig.|8]has £ « 7x 10~ 3 , which implies 
that B max = 4 mG. If the tangled field in the WCR has 
its origin in the O-star photosphere we deduce that the 
B-field at the surface of the O-star is B* ~ 100 G. This 
is comparable with the non-detection of the Zeeman ef- 
fect for the majority o f OB stars wh ich puts an upper 
limit ~ 100 G on (Mat hylll999l) . It is also worth 
noting that it is assumed that the nonthermal electron 
energy density is in equipartition with the magnetic 
and thermal energy densities. Since the energy den- 
sity of nonthermal protons, U p , is likely to be 50 x 
higher than that of the electrons, the deduced value of 
( » 7 x 10~ 3 means that U p is uncomfortably close 
to unity. If the shocks in CWBs really are this efficient 
at accelerating particles, the thermal properties of the 
shock-heated, X-ray emitting gas will be affected. It is 
interesting to note that X-ray spectra from CWBs has 
to date always been modeled and interpreted assuming 
that the shocks place an insignificant fraction of their 
energy in relativistic particles, so this finding has im- 
plications for the analysis of X-ray data. 

To appreciate how the intensity distributions from 
the models compare with the MERLIN observations 



shown in Wil liams et all (|l997), visibilities were gen- 
erated from the model intensity distributions. The re- 
sulting visibilities were then i maged and deconvolve d 
using the same procedure as in lwilliams et a l. ( 1997), 
giving the "simulated" observations shown in the upper 
panels of Fig. |9] The remarkable similarity between 
these images and the actual observations shown in the 
lower panels of Fig. |9]is striking. However, a critical 
eye reveals that the 5 GHz peak intensity of the colli- 
sion region is a little lower than that of the WR stellar 
wind in the model, and the WR star is a little too bright 
at 1.6 GHz. The WR star also appears to be a little 
more compact than the observations show, and we at- 
tribute these minor differences to our simple spherical, 
isothermal model of the stellar winds. Note that the 
OB star is not visible in the simulated observations. 
This is because the observations are essentially the in- 
tensity distribution convolved with the interferometer 
beam and the emission from the OB star is then suffi- 
ciently dispersed that it is no longer visible. 

Examination of the simulated images can help to con- 
strain the inclination angle. The synthetic observations 
shown in Fig. [9] are for an inclination of 0? At an incli- 
nation angle of 30°, the synchrotron emission is spread 
out over a much larger area and has a much lower sur- 
face brightness than in the simulation for i = 0°, and 
in the actual observations. This fact, and the spectral 
modeling, points to a fairly low value for the system 
inclination, between 0° an d 30°, though this is con- 



trary to the conclusions of Pittard et al. (2002) where 



larger inclination angles were required to recover the 
extended X-ray emission. This apparent contradiction 
may be resolved if the stellar X-ray emission is ex- 
tend ed on larger s cales than previously thought (see 



Skinner et all 2002). 



5 Additional considerations and future per- 
spectives 

In this section we briefly review processes which, 
though important, have not been discussed elsewhere 
in this paper. We r efer the interested reader to 
Folini & Waldei d2000h where some of these effects 



are discussed in greater detail, including speculations 
on their mutual interaction. Theoretical considerations 
for clumped winds and th e intriguing problem of dust 
formation are discussed in Folini & Walder (2000) and 
IWalder & Folini J2002h . 
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5.1 Radiation dynamics 

The influence of one star's radiation field on 
the dynamics of the other star's wi nd ha s been 



quantified hv Istevens & Pollockl 111994 and 



Gavl ev. Owocki & Crannied (1997, see also |Pittard, 
1998h . In close binaries the winds may be inhibited 



in their initial acceleration (since the net radiation 
flux is reduced between the stars). In WR+O binaries, 
where the WCR occurs closer to the O star than the 
WR star, the WR wind may be rapidly braked if there 
is strong coupling with the O-star radiation field. The 
opening angle of the WCR is a sensitive diagnostic 
of the braking effect, although the geometry of the 
WCR is an ambiguous constraint when the momentum 
ratio, rj, of the winds is not known. The reduction in 
pre-shock velocity and the change to the morphology 
of the WCR have obvious consequences for the X-ray 
emission, but should also impact the radio emission. 

5.2 Orbit-induced curvature of the WCR 

Orbital motion will distort the structure of the 
WCR into a large-scale spiral form, and is strik- 
ingly displayed in infrared observ ations of WR 104 
(Tuthill. Monnier & Danchi, 1999). Theoretical mod- 
els of the curva t ure range f r om th e analytical 
(Canto et al. , 1999; Tuthil l et all [2 003). to the nu- 
merical dWilliams. van der Hucht & Sooelstra. 1994), 
to the hydrodynamic al (IPittard & S tevensL [l999; 
Walder & Folinl 2003). The orbit-induced curvature 
of the WCR will affect the ob served X-ray emission 
in close binaries (IPittar d. 2002) and the radio emission 
seen from wider (especially eccentric) binaries. Theo- 
retical simulations are currently thin on the ground. 

5.3 Instabilities 

The WCR is subject to many different types of in- 
stability. These include the classical Rayleigh-Taylor, 
Kelvin-Helmholtz and Richtmyer-Meshkov instabili- 
ties which act on the surface between the two winds, 
thin-shell instabili ties which act o n isothermal shock- 
bound ed slabs ( Vishniac . 1994 : Blondin & M arks. 
1996), and the thermal instability from radiative 
cooling (e.g., Istrickland &~Blo ndin. 1995). Since 



these instabilities are not isolated from one another, 
the resulting flow may be very complex (see, e.g., 
Stevens. Blondin & PoliockL Il992h . While the degree 
of instability in theoretical models may depend on 
assumptions made in the hydrodynamical code (see 



Mvasnikov. Zhekov & Belov . 1998h . it seems certain 
that the WCR will often be dynamically unstable. This 
may lead to some X-ray dPittard & StevensL 1 1997) and 
radio variability. 

5.4 Ionization equilibrium 

It is generally assumed that the shocked gas is in 
ionization equilibrium, though in wide binaries there 
could be some departure from this (see discussion in 
Henley et al. in these proceedings). The self-consistent 
inclusion of non-equilibrium ionization has yet to be 
incorporated into hydrodynamical models of CWBs. It 
will be interesting to discover if future work in this di- 
rection is able to better reproduce the X-ray line shifts 
seen in WR 140 where inconsistencies with current 
models exist. 

5.5 Ion-electron temperature 

In wide binaries the shocks will be collisionless, and 
the immediate post-shock electron and ion tempera- 
tures may differ. This possi bility has been inve s tigate d 
in a preliminary fashion by Zhe kov & Skinned (|2000), 
where the effect on the X-ray emission from the WCR 
in WR 140 was examined. While the broad-band spec- 
tral fits to ASCA data indicated that a model with this 
effect was preferred, it is probably too soon to draw 
strong conclusions from this work and further study is 
needed. 

5.6 Thermal conduction 

There has been little work to date on thermal 
conduction in CWBs. However, the high tem- 
peratures and the steep temperature gradients in 
these systems indicate that thermal conduction 
may play an important role. Prel iminary work 
has b ee n carried out bv iMvasnikov & Zhekov 
dl998h. iMotamen. Walder & Folinil dl999h . and 
Zhe kov & Mvasnikov! d200*3l) ! where it is shown 
that efficient thermal conduction increases the X-ray 
luminosity and softens its spectrum. However, the 
presence of magnetic fields greatly reduces the thermal 
conductivity in a direction perpendicular to the field 
lines, and this has not yet been considered. This 
should be addressed in future work since we know that 
magnetic fields are present in these systems due to the 
synchrotron emission which is observed. 
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5.7 IC cooling of nonthermal particles 



While the simple models presented in Dougherty et al 
(2003) are able to reproduce both the spectrum and 
the spatial distribution of radio emission remarkably 
well, the neglect of IC cooling is a major omission. To 
realistically model the observed radio emission from 
CWBs the downstream evolution of the nonthermal 
electron energy spectrum should be calculated. This is 
being addressed in a follow-up paper (Pittard et al., in 
preparation) where, in addition, the effect of Coulom- 
bic cooling on the lower energy relativistic particles is 
also considered. 

5.8 IC cooling of thermal particles 

IC losses may be important not only for the nonther- 
mal shock-accelerated particles, but, at times, per- 
haps also for the hot p ost-shock thermal gas (see, e.g., 
White & Chen, 1995). More work is needed in this 



area to determine when this might lead to a significant 
change in the resulting thermal X-ray emission. 

6 Summary 

The wide range of parameter space covered by CWBs 
means that they are useful probes of a wide range of 
physical phenomena. X-ray observations in close bi- 
naries are a sensitive tool to study radiative braking, 
while in wide binaries X-ray emission is a useful di- 
agnostic of the electron-ion temperature equilibration 
timescale. The input of theoretical X-ray spectra into 
XSPEC via user generated "table models" has proved 
to be extremely useful in obtaining direct constraints 
of the mass-loss rates of the binary components, and 
enjoys several advantages over more traditional tech- 
niques. Such work applied to rj Carinae may in future 
determine if there is enhanced stellar mass-loss associ- 
ated with the periodic close approach of the stars. 

On the other hand, radio observations can be used to 
explore the efficiency of diffuse shock acceleration at 
densities much higher than those in other astronomi- 
cal objects with high Mach number shocks, e.g., su- 
pernovae. Preliminary work in this area indicates that 
it may be necessary to consider the effect of parti- 
cle acceleration on the shock jump conditions. An- 
other consequence of particle acceleration is that sub- 
sequent IC cooling of the relativistic particles leads 
to nontherm al emi s sion at X-ray and 7-ray en ergies, 
e.g.. lPollockl<ll987h . lBenagha & Romerol J2003h . Such 
emission may already have been detected at X-ray en- 



ergies I Rauw et al. , 2002h . and new and forthcoming 
missions like INTEGRAL and GLAST may observe this 
effect at 7-ray energies. 

Future work should address these issues and those 
identified in Sec. [5] and where possible radio and X- 
ray data should be compared against theoretical predic- 
tions from a unified model. The review of (some) im- 
portant processes and mechanisms in Section[5]demon- 
strates the richness of physical effects occurring in col- 
liding wind binaries which we are striving to under- 
stand. 
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Figure 9: Intensity distributions at 4 8 (left) and 1 .6 GHz ( right) of our WR 147 model at 0° inclination (top panels) and the 
actual observations reproduced from Willia ms et alJ (11997) (bottom panels). The simulated MERLIN observations use the 
same u-v distribution and beam sizes as the actual observations in the bottom panels. The beams are circular, of diameter 57 
and 175 mas at 4.8 and 1.6 GHz respectively. The contour levels in the simulated and actual images are also the same. The 
similarity of the simulated and actual observations is striking. 



